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Abstract. We have been proposing two evolutionary paths to Type la 
supernovae (SNe la), which are called the supersoft X-ray source (SSS) 
channel and the symbiotic channel, depending on the orbital separation 
just prior to an SN la explosion. The essential difference of our treat- 
ment is inclusion of strong, optically thick winds on mass-accreting white 
dwarfs (WDs) in the elementary processes of binary evolution when the 
mass accretion rate on to WDs exceeds a critical rate of M cr ~ 1 x 1O -6 M0 
yr _1 . Once optically thick winds begin to blow from the WD, the binary 
can avoid forming a common envelope in some cases. We call this accre- 
tion wind. So that the WDs are able to grow up to the Chandrasekhar 
mass and explode as an SN la, showing SSS or recurrent nova phenomena 
in the way to SNe la. Thus, the accretion wind process of WDs can open 
new channels to SNe la. We have modeled the LMC supersoft source RX 
J0513. 9-6951 as an example of the systems in the accretion wind phase. 
Further inclusions of two other elementary processes make the channels 
much wider; these are the case BB mass transfer in the SSS channel and 
the strong orbital shrinkage during the superwind phase of the primary 
star in the symbiotic channel. As a result, the estimated birth rate of 
SNe la via these two channels becomes compatible with the observation 
in our Galaxy. Interestingly, the U Sco and T CrB subclasses of recur- 
rent novae can be naturally understood as a part of evolutionary stages 
in these two SSS and symbiotic channels to SNe la, respectively. Thus 
we have a unified picture of binary evolutions to SNe la, luminous SSS, 
and recurrent novae. 



1. Introduction 

Type la supernovae (SNe la) are one of the most luminous explosive events of 
stars. Recently, SNe la have been used as good distance indicators that provide 
a promising tool for determining cosmological parameters (Perlmutter et al. 
1999; Riess et al. 1998) because of their almost uniform maximum luminosities. 
However, one of the most important unresolved problems on SNe la is that we 
do not know the exact progenitor systems of SNe la. 

It is widely accepted that the exploding star itself is an accreting white dwarf 
(WD) in a binary (Nomoto 1982; Nomoto, Thielemann, & Yokoi 1984). However, 
the companion star (and thus the observed binary system) is not known. Several 
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Figure 1. Left: Optically thick winds blow from mass-accreting 
white dwarfs when the mass transfer rate exceeds a critical rate, i.e., 
Macc > M cr . Here, we assume that the white dwarf accretes mass 
from equatorial region and, at the same time, blows winds from polar 
regions. Right: Response of white dwarfs to mass accretion is illus- 
trated in the white dwarf mass and the mass accretion rate plane, i.e., in 
the MwD-M acc plane. Strong optically thick winds blow above the line 
of M acc > M cr . There is no steady state burning below M acc < M st( j. 
Instead, intermittent shell flashes occur. The envelope mass at which 
hydrogen ignites to flash is also shown (taken from Fig. 9 of Nomoto 
1982). 



objects have ever been considered, which include merging double white dwarfs 
(Iben & Tutukov 1984; Webbink 1984), recurrent novae (Starrfield, Sparks, & 
Truran, 1985), symbiotic stars (Munari & Renzini 1992) etc (Livio 2000, for 
recent summary). 

Among these models, the systems of double white dwarfs (the double de- 
generate model) are not theoretically supported; when two WDs merge, carbon 
ignites not at the center but near the surface layer of the C+O WD; so that 
the C+O WD is converted into an O+Ne+Mg WD (Saio & Nomoto 1998). If 
the total mass of the merger exceeds the Chandrasekhar mass, the merged WD 
does not explode as an SN la but collapses to a neutron star (Nomoto & Kondo 
1991). 

On the other hand, the single degenerate systems has become more promis- 
ing, in which a C+O WD gains mass from a non-degenerate companion star 
and grows its mass up to the Chandrasekhar mass. Such promising progenitors 
include luminous supersoft X-ray sources (SSS), recurrent novae (RNe), and 
symbiotic stars. However, the most important problem in these single degener- 
ate systems is that we did not know the mechanism how to grow the mass of the 
WDs. Instead, there are effective mass-losing mechanisms such as nova ejection 
and common envelope evolutions that get rid of the mass accumulated on the 
WD. Here, I review and summarize our recent developments of the mechanisms 
that enable to grow the mass of WDs in the single degenerate systems and give 
a unified picture of evolutions of SSSs, recurrent novae, and SNe la. 
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2. Accretion wind evolution 

First, we briefly explain how the mass of WDs in the single degenerate systems 
grow toward the Chandrasekhar mass limit. The reason that the standard binary 
evolution theory fails to explain the evolutionary path to SNe la and recurrent 
novae lies mainly in the theoretical prediction of a common envelope formation 
and the ensuing spiral-in process, because these processes inhibit the growth 
of white dwarfs in binary systems. As discussed in many previous papers on 
binary evolution, it has been widely accepted that a hydrogen-rich envelope on 
a mass-accreting white dwarf expands to a red giant size when its mass accretion 
rate, M acc , exceeds a critical limit, i.e., 

M cr = 8.3 x 10" 7 (^ - 0.4) M Q yr-\ (1) 

(see, e.g., Fig. 9 of Nomoto 1982) and easily forms a common envelope. Once a 
common envelope is formed, two stars begin to spiral-in each other as a result 
of viscous drag, thus producing a double degenerate system (Iben & Tutukov 
1984; Webbink 1984). 

However, Hachisu, Kato, & Nomoto (1996) found that white dwarfs begin 
to blow strong winds (v w ind ~ WOO km s _1 and M w j n d > 1 x 1CT 6 M yr _1 ) when 
M acc > M cr , as illustrated in Fig. 1, thus preventing the binary from collapsing. 
The binary does not spiral-in but keeps its separation. Hydrogen burns steadily, 
and therefore the helium layer of the white dwarf can grow at a rate of Mn e ~ 
M cr . The other transferred matter is blown in the wind (M w j n d ~ M acc — M CI ), 
where M w - mc \ is the wind mass loss rate. We call this accretion wind because it 
begins to blow when the accretion rate exceeds the critical rate. 

Therefore, we have to revise Fig. 9 of Nomoto (1982), which has been 
used widely in binary evolution scenarios to show the basic processes of mass- 
accreting white dwarfs. The revision to Nomoto's Fig. 9 is shown in Fig. 1. The 
most important difference of Fig. 1 from Nomoto's Fig. 9 is that white dwarfs 
increase their masses in a much wider parameter region. In our new diagram, 
the status of a mass-accreting white dwarf is specified by the following three 
phases: (1) accretion wind phase (M acc > M cr ); (2) steady shell-burning phase 
(M st d < M acc < Mcr)', and (3) intermittent shell-flash phase (M acc < M st d), 
where M st( j means the lowest limit of the mass accretion rate for steady hydrogen 
shell-burning. The new growing region of white dwarfs, M acc > M st d, is much 
wider than old Nomoto's (1982) narrow strip, M st d < M acc < M cr . 

2.1. RX J0513. 9-6951 as an example of accretion wind phase 

As an example of the accretion wind phase, we have modeled long-term optical 
light curves of the LMC supersoft X-ray source RX J0513. 9-6951 based on our 
wind model of mass-accreting white dwarfs. The source shows a quasi-periodic 
variability, the period of which is ~ 170 d and the amplitude of which is ~ 1 
mag (Alcock et al. 1996). The source has a relatively longer (~ 130 d) optical 
high-state (V ~ 16.2) and a shorter (~ 40 d) optical low-state (V ~ 17.0) with 
the transition time of several days. X-rays were observed only in the optical low- 
states. This indicates that a strong wind blows only in the optical high-state 
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Figure 2. Left: Geometrical configuration of our RX J0513. 9-6951 
model in the optical high state (wind phase). The cool component is a 
slightly evolved MS companion filling up its inner critical Roche lobe. 
The hot component, WD, blows a strong wind. The accretion disk is 
blown off outward and has a larger size than the no- wind phase. Right: 
Geometrical configuration in the optical low state (no- wind phase). 
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Figure 3. Light curve is plotted against time for a 1.35M Q white 
dwarf model together with the MACHO data (Alcock et al. 1996). We 
assume a larger size of the accretion disk (1.0 times the Roche lobe 
size) in the optical high state (the wind phase) but a smaller size (0.8 
times the Roche lobe size) in the optical low state (the static, no-wind 
phase). The transition from the optical high to low state corresponds 
to the cease point of the strong wind (large open circles). 
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because X-rays are self-absorbed by the wind. To examine the idea of the strong 
wind, we have calculated envelope solutions of mass-accreting WDs and found 
that optically thick winds really occur when the mass accretion rate exceeds a 
critical rate of M cr for J%d > O.7M even in a low metallicity of Z = 0.004 of 
LMC. 

We attribute the 100 — 150 d bright phase to our strong wind phase with an 
expanded white dwarf photosphere, whereas the ~ 40 d less bright phase to our 
static (no-wind) phase with a shrinking envelope of the WD. The duration and 
the brightness of the optical high state is explained by the model of a WD mass, 
M^d = 1-35 Mq, with an extending accretion disk edge over the Roche lobe 
size during the wind phase as illustrated in Fig. 2. Irradiation of the accretion 
disk by the WD photosphere plays an essential role to reproduce the optical 
light curve because a large part of the optical light comes from the accretion 
disk. 

The switching mechanism from the optical high to low state is understood 
as follows: Strong winds from the WD themselves make the separation increase, 
because the high velocity wind do not carry away large angular momentum. The 
winds further strip off the very surface layer of the companion star. As a result, 
the mass transfer rate decreases and it finally stops the strong winds. After the 
wind stops, the surface of the companion recovers in a thermal timescale of the 
envelope and the mass transfer rate increases up and strong winds begin to blow 
again. Thus, the semi-periodic cycles repeat. 

Light curve fittings indicate that the optical high state is initiated by a 
sudden increase in the mass transfer rate of M a cc ~ 1 x 10~ 5 Mq. Then, the 
envelope of the white dwarf expands to R p h ~ 0.1 R G , not exceeding over the 
accretion disk size of ~ 2 Rq. Our light curve fittings are shown in Fig. 3. In 
our models, about 40% of the transferred matter can be accumulated on the 
white dwarf but the rest is blown off in the wind. This suggests that it takes 
~ 1 x 10 5 yr for the white dwarf of 1.35 M to grow to 1.38 M Q and to explode 
as a Type la supernova. 

2.2. Two evolutionary paths to RNe and SNe la 

Based on the new mechanism of the accretion wind evolution as shown in Fig. 
1, Hachisu, Kato, & Nomoto (1999a) and Hachisu et al. (1999b) have proposed 
two paths to SNe la (see also, Li & van den Heuvel 1997). In these two paths, 
WDs can accrete mass continuously from the companion and grow at a rate of 
M cr of eq. (1). Therefore, candidates for SN la progenitors are systems at the 
final stages on these paths; i.e., one is a supersoft X-ray source (SSS) consisting 
of a WD and a lobe-filling, more massive, main-sequence (MS) companion as 
shown in the late stage of Fig. 4 (the WD+MS systems; see also Li & van 
den Heuvel 1997) and the other is a symbiotic star consisting of a WD and a 
mass-losing red giant (RG) as shown in the late stage of Fig. 5 (the WD+RG 
systems, see also Hachisu et al. 1996). Both the systems contribute to main 
parts of the SN la birth rate (Hachisu et al. 1999a, 1999b). 

Hachisu et al. (1999a, 1999b) have followed many binary evolutions until 
SN la explosions. The final outcome is summarized in Fig. 6 for the initial 
orbital period, Pq (in units of days), and the initial donor mass, M^o, which is 
representing Mrg,o (the initial mass of the red-giant component), or Mms,o (the 
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Figure 4. SSS channel. Left: Early evolutionary path through the 
common envelope evolution to the helium (Case BB) mass transfer. 
Right: Late evolutionary path to an SN la explosion in our wind 
model. 



initial mass of the slightly evolved main-sequence component). In this figure, 
the initial mass of the white dwarf is assumed to be M\vd,o = 1-OM . The 
right-hand region, longer orbital periods of Pq ~ 100 — 1000 d, represents our 
results for the WD+RG system. In the left-hand region, shorter orbital periods 
of Pq ~ 0.3 — 10 d, are shown the results for the WD+MS system. 

We plot the initial parameter regions of binaries (thin solid curves in Fig. 
6) that finally produce an SN la, which are corresponding to stage (F) in Fig. 
4 or stage (D) in Fig.5. Left diagram in Fig. 6 shows such two regions of the 
WD+MS systems (SSS channel) and the WD+RG systems (symbiotic channel). 
We have added in this figure the final parameter regions (thick solid curves) 
just before an SN la explosion occurs. Starting from the initial region encircled 
by the thin solid curves, binary systems evolve and can explode as an SN la in 
the regions enclosed by the thick solid curves. Hachisu & Kato (2001) showed 
that some recurrent novae lie on the boarder of their parameter region that can 
produce an SN la (see also Kato 2001 in this volume). Among six recurrent novae 
with known orbital periods, five fall in the final regions of SN la progenitors just 
before the explosion. Only T Pyx is far out of the regions of SNe la. 



3. Case BB mass transfer 

It has been shown that the ejecta of the recurrent nova U Sco are helium-rich 
as much as He/H~ 2 by number (Williams et al. 1981). The accretion disk of 
U Sco is also suggested to be helium-rich (Hanes 1985). Many of the SSSs have 
a stronger He II line than Balmer lines, thus also suggesting that the companion 
of the WD+MS systems (SSS channel) is helium-rich. In this section, therefore, 
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Figure 5. Symbiotic channel. Left: Early evolutionary path 
through the first common envelope evolution to a symbiotic system. 
Right: Late evolutionary path through the symbiotic channel to an 
SN la. 



we describe the reason that the cool component of the WD+MS systems has a 
helium-rich envelope. 

Previous works consider only the initial systems consisting of a more massive 
AGB star with a C+O core and a less massive main-sequence star. This system 
undergoes a common envelope evolution and it finally yields a binary system 
of a mass-accreting C+O WD and a lobe-filling MS/subgiant star. Thus, they 
do not include another important evolutionary path, in which a more massive 
component fills up its inner critical Roche lobe when it develops a helium core 
of ~ 0.8 — 2.0M Q in its red giant phase as shown in the early stage of Fig. 4. 

We consider, for example, a close binary with the primary mass of Mi j = 
7Mq and the secondary mass of M<i^ = 2M with the initial separation of 
cij ~ 80 — 600 i? Q . When the primary has evolved to a red giant developing a 
helium core of M^He ~ 1.0— 1.6M , it fills up its inner critical Roche lobe. Then, 
the binary undergoes a common envelope (CE) evolution, which yields a much 
more compact close binary of a naked helium star of M^hc ~ 1.0 — 1.6M Q and a 
main-sequence star of M2 = 2M & with the final separation of a+- 5 cE ~ 4 — 50 R Q . 
Here, we assume a relation of 



Of,CE 



«CE 



Mi ;C 



Mo 



Mi i - Mi, 



(2) 



with the efficiency factor of «ce = 1-0 for the common envelope evolutions. 

After the common envelope evolution, the radii of the inner critical Roche 
lobes of the primary and the secondary become Rl ~ O^a^cE and R\ ~ 
O^a^cE, respectively. Since the secondary radius should be smaller than its 
inner critical Roche lobe, i.e., 1.7-R0 < R<i < R\ (1.7 Rq is the radius of a 2M 
ZAMS star), the initial separation should exceeds 80 R & . The upper bound 
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of the initial separation is obtained from the maximum radius of the 7M star 
which has formed a helium core, i.e., a« ps 2i?i >max < 2 x 300 Rq. Thus, 
the allowable range of the initial separations is 80 Rq < a« < 600 Rq for a 
pair of M M = 7M and M 2)i = 2M . Then, we have ~ 1.5 - 18 R Q and 
i?2 ~ 1.7 — 20 i?© after the first common envelope evolution. 

After the hydrogen-rich envelope is stripped away and hydrogen shell burn- 
ing vanishes, the naked helium core contracts to ignite a central helium burning 
and becomes a helium main sequence star. When the helium star forms a C+O 
core, its helium envelope expands to fills its inner critical Roche lobe again. The 
helium is transferred stably to the secondary on an evolutionary time scale of 
tev ~ 10 5 yr because the mass ratio is smaller than 0.79 (q = M1/M2 < 0.79). 
The resultant mass transfer rate is Mi ~ 1O~ 5 M yr" 1 . No common envelope 
is formed for such a low rate. After the helium envelope is lost, the primary 
becomes a C+O WD of M +o ~ 0.9 — 1.1M and the separation increases 
by 10%-40%, i.e., a f)He ~ (1-1 - L4)a f)C E ~ (4 - 7O)i? . Here, we assume 
the conservations of the total mass and angular momentum during the helium 
mass transfer to obtain the separation after the helium mass transfer, a^He- 
The secondary receives almost pure helium matter of AMn e ~ 0.1 — O.5M to 
form a helium-enriched envelope. The secondary's hydrogen content decreases 
to X ~ 0.6 if helium is completely mixed into the star. For the 9 + 2.5M case, 
much more helium is transferred, but much less helium for the 6 + 2M case. 

4. Shrink of orbital separation during superwind phase 

Based on the population synthesis analysis, Yungelson & Livio (1998) claimed 
that almost no realization frequency is derived for the original WD+RG model 
(Hachisu et al. 1996). In this section, we point out that very wide binaries with 
initial separations of a\ > 1500 Rq, which were not included in their analysis, 
are essentially important in our SN la modeling. 

The more massive component (the initial mass of M\^) of a binary first 
evolves to a red giant (AGB stage) and fills its inner critical Roche lobe. After a 
common envelope phase, the more massive component leaves a C+O WD, and 
the separation of the binary decreases by a factor of aj/ai ~ 1/40 — 1/50 from 
eq.(2) for M WD ~ 1 M Q and M 2 ~ 1 M because a ~ 1 M Q WD descends 
from a main-sequence star of ~ 7 — 8 M . Yungelson & Livio (1998) 
assumed that the separation of interacting binaries is smaller than ai < 1500 Rq. 
Then, the widest binaries have separations of aj < 30 — 40 Rq after common 
envelope evolution. The orbital period is Pq < 20 days for Mwd,o ~ 1 M Q 
and Mrg,o ~ 1 Mq. There is no SN la region of the WD+RG systems for 
Pq < 20 days as seen from Fig. 6. Thus, Yungelson & Livio concluded that we 
cannot expect any SN la explosions from the right-hand SN la region (WD+RG 
system) in Fig. 6. 

The reason that Yungelson & Livio's modeling fails to reproduce our SN 
la progenitors is their assumption of a\ < 1500 Rq. In what follows, we show 
that WD+RG binaries wide enough to have Pq ~ 100 — 1000 d are born from 
initially very wide binaries with a, ~ 1500 — 40000 Rq. 

A star with a zero-age main-sequence mass of < 8 M ends its life by 
ejecting its envelope in a wind of relatively slow velocities (v ~ 10 — 40 km s" 1 ) 
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Figure 6. Left: Regions producing SNe la are plotted in the log P— 
Md (orbital period — donor mass) plane for the WD+MS system (left) 
and the WD+RG system (right). Here we assume an initial white 
dwarf mass of 1.0M Q . Currently known positions of each recurrent 
nova are indicated by a star mark (*) for U Sco (Hachisu et al. 2000), 
a filled triangle for T CrB (Belczyhski Sz Mikolajewska 1998), an open 
rectangle for T Pyx (Patterson et al. 1998), but by arrows for the 
other three recurrent novae, V394 CrA, CI Aql, and RS Oph, because 
the mass of the companion is not yet available explicitly. Right: The 
accretion winds of WDs open new two channels to SNe la while the 
two evolutionary processes of the case BB mass transfer and the orbital 
shrinkage during the superwind phase enlarge the initial phase volume 
of SN la progenitors for the SSS channel and the symbiotic channel, 
respectively. 

as shown in Fig. 5. These wind velocities are as low as the orbital velocities of 
binaries with separations of aj ~ 1500 — 40000 R & for a pair of ~ 7 M Q 
and M2 t i ~ 1 Mq. When the wind velocity is as low as the orbital velocity, we 
have the separation shrinking due to effective loss of the angular momentum by 
the outflowing gas, i.e., a/a ss 2M1/M2 < (see Hachisu et al. 1999a). 

Once the binary system begins to shrink, its evolution becomes similar to 
a common envelope evolution (see stages B-D in Fig. 5). Thus, the separation 
is reduced by a factor of 1/40 - 1/50, i.e., a f ~ 30 - 1000 R Q for M M ~ 7 M Q 
and Mi,% ~ 1 Mq. The orbital period becomes Pq ~ 15 — 3000 days for a pair 
of Mwd,o ~ 1 Mq and M2 ~ 1 M & (stage D in Fig. 5). These initial sets of 
the parameters are very consistent with the initial conditions of our WD+RG 
progenitor systems. 

5. Type la supernova rates 

We estimate the SN la rate from the WD+RG and WD+MS systems in our 
Galaxy by 



where q, A, and M are the mass ratio (q = MijM\ < 1), the initial separation 
in units of solar radius, and the primary initial mass in solar mass units (from 




-1 



(3) 
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cq. [1] of Iben k Tutukov 1984). The SN la rate is essentially the same as the 
volume of the initial phase space consisting of q, A, and M as shown in Fig. 6. 

The accretion winds of WDs open new two channels to SNe la. The case BB 
mass transfer in the SSS channel extends the phase volume towards the low mass 
ratio region as shown in Fig. 6. On the other hand, the orbital shrinkage during 
the superwind phase in the symbiotic channel enlarges the phase space towards 
the large initial separation also shown in Fig. 6. Therefore, the estimated rate of 
WD+RG/WD+MS systems becomes close to the observed rate in our Galaxy, 
i.e., v ~ 0.003 yr -1 , which comes from the SSS channel of i^ms ~ 0.001 yr _1 and 
the symbiotic channel of ^rg ~ 0.002 yr _1 . 
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